0++ Stat es in a large-Nc Regge approach 
Enrique Ruiz Arriola 

Departamento de Fi'sica Atomica, Molecular y Nuclear 
and Institute Carlos I de Fisica Teorica y Computacional, 
Universidad de Granada, E-18071 Granada, Spain 

AND 

WOJCIECH BrONIOWSKI 

The H. Niewodniczaiiski Institute of Nuclear Physics PAN, PL-31342 Krakow, Poland, and 
Institute of Physics, Jan Kochanowski University, PL-25406 Kielce, Poland 

(Received 12 October 2011) 

Scalar-isoscalar states (J^*^ = O"*""*") are discussed within the large-iVc Regge approach. We 
find that the lightest /o(600) scalar-isoscalar state fits very well into the pattern of the radial 
Regge trajectory where the resonance nature of the states is advantageously used. We confirm 
the obtained mass values from an analysis of the pion and nucleon spin-0 gravitational form 
factors, recently measured on the lattice. We provide arguments suggesting an alternating meson- 
glueball pattern of the 0++ states, which is supported by the pseudoscalar-isovector 0'' excited 
spectrum and asymptotic chiral symmetry. Finally, matching to the OPE requires a fine-tuned 
mass condition of the vanishing dimension-2 condensate in the Regge approach with infinitely 
many scalar-isoscalar states. 



PACS numbers: 12.38.Lg, 11.30, 12.38.-t 



1. Introduction 

The goal of this talk is to discuss various intriguing aspects of the spectrum of scalar-isoscalar 
states. Approaches developed in recent years may shed new light on this long-elaborated problem 
in hadronic physics. 

Subsequent hadron resonances listed in the Particle Data Group (PDG) tables increase their 
mass up to the upper experimental limit of 2.5 GeV, while their width remains bound within 
500 MeV. In Fig. [T] we show separately the widths of all baryons and mesons listed in the PDG 
tables \1\ as functions of the mass of the state. One naturally expects that broad resonances, 
i.e., with r ~ m, escape phenomenological analysis; even if they existed, they might be missing 
from the PDG as difficult to assess experimentally. Note, however, that with the exception of 
the notorious /o(600) resonance and a few baryon and meson states, the ratio is bound by the 
line r/m ~ 1/3 (the dashed line in Fig. [2]) . 

A natural and model-independent framework to understand this feature is provided by the 
limit of large number of colors in QCD. Indeed, in this large- iVc limit, with g^Nc fixed, baryons 
are heavy with mass m = 0{Nc) and width T = 0{N^) [2|[3], while mesons and glueballs are 
stable, with mass independent of Nc, namely m = 0{N^), and width T suppressed as 1/Nc and 
respectively. This means that P/m is suppressed (see, e.g., ji] for a review). In particular, 
one has F/m ~ for mesons and baryons, while P/m ~ for glueballs 

* Presented by ERA at Miniworkshop Understanding Hadronic Spectra, Bled (Slovenia) 3-10 July 2011 
^ Fig. [l] suggests that it is reasonable to assume that excited states in the spectrum follow a more accurate 
large- A'^c pattern than the ground state. 
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Fig. 1. Width of all baryons (left panel) and mesons (right panel) listed in the PDG tables [l] as a function of 
the hadron mass (in GeV). The surface of each point is proportional to the (2 J + 1) spin degeneracy, while the 
intensity is proportional to the isospin degeneracy {21 + 1). 



all baryons 



0.0 
0.0 



1.0 1.5 

m [GeV] 



all mesons 



1 

0.5 



1.0 1.5 

m [GeV] 



Fig. 2. Same as Fig. [l] for the width/mass ratio. The large-A'^c limit predicts F/m ~ 0{Nc ^). The dashed 
horizontal line corresponds to 1/3. The gray bands reflect the uncertainties of the fit F/m = 0.12(8). 



On the other hand, at high excitations mesons are expected to resemble strings of length 
/ ~ m/(T, with a denoting the string tension (see, e.g., Ref. [s] and references therein). Actually, 
the decay rate of a string per unit time, F, is estimated to be proportional to the length (6|[7], 
which immediately yields constant F/m |^ In Fig. [2] we show the ratio of widths to masses of all 
baryons and mesons listed in the PDG tables ll], plotted as functions of the mass of the state. 
If we compute the average weighted with the (2 J + 1) spin degeneracy and its spread, we find 
(a are the remaining quantum numbers, including isospin) 

^-E(2^+ 1)^ = 0-12(8) (1) 

iiL I II' J, a 

for both baryons and mesons! This rather small ratio, complying to the large-A'^c and string- 
model arguments, suggests that this is a generic feature of the hadronic spectrum rather than a 
lack of experimental ability to resolve too broad states. One may thus assert in this regard that 
the PDG spectrum is fairly complete. 

2. Masses of resonances 

A resonance may be interpreted as a superposition of states with a given mass distribution, 
approximately spanning the m it F/2 interval. Of course, the shape of the distribution depends 
on the particular process where the resonance is produced, and thus on the background. The 
rigorous quantum-mechanical definition of the resonance corresponds to a pole in the second 



^ The argument directly carries over to baryons, treated as a quark-diquark string. 



Ruiz-Arriola printed on October l^, 2011 



3 




Fig. 3. Radial Regge trajectory corresponding to the squared mass of all = 0^"'' scalar-isoscalar and = 
^ pseudoscalar-isovector states listed in the PDG tables |T|. The four heaviest 0^^ and two ^ states are 
not yet well established and are omitted from the PDG summary tables. The error bars correspond to the errors 
in the determination of the square of mass as Am^ = mF with F from [l]. The straight line is the result of our 
joined fit. Labels of 0"'"'' states are above their mark whereas labels of 0"' states are below their mark. 



Riemann sheet in the partial- wave amphtude of the corresponding decay channel, which becomes 
independent of the background. However, although quoting the pole is highly desirable, with a 
few exceptions this is not what one typically finds in the PDG. 

As a matter of fact, several definitions are used: pole in the second Riemann sheet, pole 
in the iC-matrix, Breit-Wigner resonance, maximum in the speed plot, time delay, etc. (see, 
e.g., (8||9]). Clearly, while all these definitions converge for narrow resonances, even for broad 
states we expect the masses to be compatible within their corresponding m it r/2 intervals. As 
mentioned, the values listed in the PDG for a given resonance correspond to different choices 
and/or processes, but mostly the results are compatible within the estimated width differences. 

The lowest resonance in QCD is the O"*""*" state /o(600) or the cr— meson. It appears as a 
complex pole in the second Riemann sheet of the vrvr scattering amplitude at ^/s^ = mu — 
with ma = 44llf and T„ = 5Utll MeV (lo| (see also Ref. (111). While these are 
remarkably accurate, it is unclear whether these numbers can be directly used in hadronic 
physics. An analysis of the role played by the u as a correlated 27r exchange in the central 
component to the NN force shows that the complex-pole exchange does not accurately describe 
this effect in the range 1 fm < r < 5 fm, but prefers a value in between the pole and the Breit- 



Wigner approximation, rua = 600(50) MeV 12 . The Breit-Wigner parameters are (mo-. To-) 



(841(5), 820(20)) MeV (l3j. For the time delay method we obtain {ma,Ta) ~ (475,630) MeV. 
As we see, the different determinations agree within the wide mo- it ro-/2 interval. 



3. Scalar Regge spectrum 

Higher 0^+ states listed in the PDG [l] follow roughly the general pattern of increasing mass 
but not their width. Radial and rotational Regge trajectories were analyzed in Ref. 14 . For 



16 



trajectories could then be identified, including three states 
T?! (see also [l8| ) we have analyzed all the 0"*"^ states which 
appear in the PDG tables (see Fig. 3^and found that all fit into a single radial Regge trajectory 
of the form 



scalar states [15] two parallel radial 
per trajectory. In a recent work 



Msinf 



-n + m„ 



(2) 



The mass of the a state can be deduced from this trajectory as the mass of the lowest state. 
The resonance nature of these states suggests to use the corresponding half- width as the mass 
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uncertainty in the fit 

Minimization yields x^/DOF = 0.12, with 

a = 1.31(12) GeV^ m„ = 556(127) MeV. (4) 

Formula ^ is actually equivalent to two parallel radial Regge trajectories with the standard 
slope, 

Ms,^{nf = an + ml, (5) 
Ms,+ {nf = an + ml + (6) 

where a = 2TTa, and a is the string tension associated to the potential V{r) = ar between 
heavy colored sources. The value ^/(T = 456(21) MeV obtained from our fit agrees well with 
lattice determinations of ^/a = 420 MeV [l9]. Of course, one expects some of these states to 
correspond eventually to glueballs. However, there seems to be no obvious difference between 
mesons and glueballs, as far as the radial Regge spectrum is concerned. Note that Casimir 
scaling suggests that the string tension is fTgiuebaii = fcmeson, but this holds in the case of fixed 
and heavy sources. The fact that we have light quarks might explain why we cannot allocate 
easily the Casimir scaling pattern in the light-quark scalar-isoscalar spectrum. 

4. Interpolating fields 

For scalar states a measure of the spectrum is given in terms of the (gauge and renorm 



invariant) trace of the energy momentum tensor 20 



= e = ^G^^^G% + [1 + jrnia)] qq. (7) 

q 

Here /3(a) = fi'^da/d^'^ denotes the beta function, a = j (47r) is the running coupling constant, 
7m (a) = dlogm/dlog/x^ is the anomalous dimension of the current quark mass m^, and GJJj^ is 
the field strength tensor of the gluon field. 

It is interesting to consider the situation of massless quarks, where things become somewhat 
simpler. Then, we have in principle two scalar operators with smallest canonical dimensions, the 
gluon and the quark qq. While these two operators are both scalars, they are chirally even 
and odd, respectively, i.e., under the q — t- 755 transformation. Because the chiral symmetry is 
spontaneously broken, there is some mixing between G^ and qq. These operators connect scalar 
states to the vacuum through the matrix element 

(0|6|n) = rr?Jn , (Olg^ln) = m„c„ . (8) 

The two-point correlators read 

ne©((?) = i [ dW^-(o|r{e(x)G(o)}|o) = V-|^ + c.t., (9) 

/ ^-^ m„ — q^ 

n " 

nes('?) = i I d''xe'^-^Q\T{e{x)qqm |0) = + e.t. , (10) 

J „ q 

I^ssiq) = i I d^xe*^-"(0|r{g-g(x)g-g(0)}|0) = j;^^ + c.t., (11) 
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where in the r.h.s. we saturate with scalar states and c.t. stands for subtraction constants which 
can chosen as to replace — )• in the numerator. In that scheme, in the large —q^ ^ Aqcd 



limit, a comparison with the Operator Product Expansion (OPE) [21-23 leads to the matching 
conditions 

neG(g') = g^Coiog(V) + -.., 

nss(g') = q^C'^\og{-q^) + ..., 

nes(g') = C^,'(gg)log(V) + ..., (12) 

where Co = -(2/3(Q)/a7r)2, Cq = -3/(87r^) and Cq = -2j3{a)/a'K. As we see, qq and do 
not mix at high q^ values, a consequence of asymptotic chiral symmetry [20^ . In these limits the 
sums over n can be replaced by integrals, whence the following asymptotic conditions are found: 

frJ (dml/dn) Cq , c^/ {dml/dn) C'q , Cnfnml/ {dm^Jdn) -f C'i^{qq) . (13) 

We see that the first two conditions are incompatible with the third one if increases for large 
n, as is the case of the data. However, if we group the states in two families, as suggested by 
Eqs. ^ and Fig. [sj we get a compatible solution 

Cn-,fn,+ const c„,+ , /„ _ const/m^. (14) 

This is equivalent to assuming an asymptotically alternating pattern of mesons and glueballs, 
coupling to chirally odd and even operators, qq and G^, respectively. Since asymptotically 

~ an/2, we find c„^_/c„_+ and fn,+/fn - ~ ns. Of course, this is not the only solution. 

The situation described above suggests the existence of a hidden symmetry in the O"*""*" sector. 
In our case we could think of the 75-parity (which becomes a good quantum number for excited 
states) as the relevant symmetry which makes the difference between the chirally even and odd 
states. This, however, only explains the fact that asymptotically the slopes of the + and — 
branches are the same, but not why the intercepts accurately differ by half the slope. 



5. The holographic connection 

To further elaborate on this intriguing point of the accidental degeneracy, let us consider the 
one-dimensional harmonic oscillator with frequency u, as an example; all states ipn{z) with the 
energy En = hoj{n + l/2) can be separated into parity even and parity odd states, satisfying the 
conditions '4'n,±{z) = ip2n{z) and '4^±,n{—z) = ±^-i-_„(z), respectively, and having the energies 
Ej^^n = 2huj{n + 1/4) and E-^n = 2hu){n + 3/4). These formulas display twice the slope of En- 
Thus, given the states with energies -E+.n and we might infer that parity was a hidden 

symmetry of a Hamiltonian explaining the correlation between the slope and intercepts. 

In the relativistic case the argument can also be made in a suggestive manner. Let us 
consider the Klein-Gordon action for infinitely many bosons in four dimensions, described with 
fields (pnix) of masses 

S=\f d^xY, [d^(t>nd^.(t>n - mlcfi] . (15) 

We assume the spectrum of the form = an + m^. Next, we can introduce the five-dimensional 
field 0(x, z) = Yin 4>n{x)'4jn[z), with ipniz) fulfilling the auxiliary Sturm-Liouville problem in the 
variable < z < 00, 

- dz [p{z)dztpn{z)] + q{z)7pn{z) = mlp{z)ipniz) , (16) 

where the functions are orthogonal with respect to the weight function p{z), provided suitable 
boundary conditions 



p{z) {'lpn{z)'>pm{z) - ■tpn{z)^p'„i{z)) \z=0 = 



(17) 



6 



Ruiz-Arriola printed on October 14, 2011 



and ^n(oo) — t- are fulfilled. The action can then be written as 

S=]^l (fx dz [p{z)d^'<t)d^(^ - p{z){dA? - q{z)(t>^] (18) 







after some integration by parts in the variable z. This action can be written as a five-dimensional 
action with a non-trivial metric |^4], featuring the AdS/CFT (soft-wall) approach (see f25] and 
references therein), with the extra dimension z playing the role of a holographic variable and 
the orthogonal set of functions ipn{z) denoting the corresponding Kaluza-Klein modes. Clearly, 
z has the dimension of length, suggesting that z — )• corresponds to the ultraviolet and z ^ oo 
to the infrared regime. 

Turning to Eq. ([16]) , we may take the standard Harmonic oscillator Schrodinger-like equation 

{p{z) = p{z) = 1, q{z) = U{z) = a^zVlG) 

- C(^) + T^a^^Vn(2:) = mliljniz) (19) 
16 

and obtain for the regular solutions at infinity the result 

r/'3 



where T{x) is the Euler Gamma function, which is meromorphic and have simple poles at 
X = 0, —1, —2, . . .. The solutions fulfilling the Dirichlet, ipniS^) = 0, and Neumann, V'n(O) = 0' 
boundary conditions, respectively, have the masses 

mi .^ = an + I , = an + ^ , (21) 

which can be merged into one single formula 

ml = -^{2n+l). (22) 

This yields rria = mjy/-v/3 = 566 MeV and, for the string tension, a = m y^2/37r = 450 MeV 
with rufg = 980 MeV, quite reasonable values. 

In this approach the symmetry in the scalar spectrum corresponds to a parity symmetry 
in the holographic z variable ^n{—z) = =b'0„(z). Note that usually the holographic variable z 
is taken to be positiv^ but if we extend it to —oo < z < oo, we may define a holographic 
superfield containing two different and orthogonal modes. Otherwise, in the interval < z < oo 
the Dirichlet and Neumann modes are not orthogonal to each other. 



6. Pseudoscalar mesons and chiral symmetry 

Discerning the nature of the a state has been a recurrent pastime for many years. As is well 
known, glueballs are more weakly coupled to mesons, 0{1/Nc), than other mesons, 0{l/\/Nc)- 
The minimum number of states, allowed by certain sum rules and low energy theorems, is just 
two. In Ref. [l6] we undertake such an analysis which suggests that /o(600) (denoted as a) is 
a qq meson, while /o(980) (denoted as /o) is a glueball. This is supported by the rather small 



^ This is supported by the hght-front interpretation of Brodsky and de Teramond 26 , where the holographic 
variable is the polar coordinate of a two dimensional vector, z = |("| and = b^x(l — x), with fo denoting the 
impact parameter and x the momentum fraction of the quark. This interpretation yields a two dimensional 
potential U{Q = + 2k?(L + S — 1) with J — L + S which, when passing to the polar variable z, generates 
the usual centrifugal term {L^ — l/4)/z^ not present in our discussion, yielding M^ j^ g = Ak^Iu + L + 5/2) 



which for J — and L — I resembles Eq. ( 22 1 
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width ratio, which yields ^f/Ta ~ (g'j^^m'j) / (g'^^^m^) ~ 1/Nc, thus for rria ~ 0.8 MeV ^ ruf 

the ratio ga-w-K / Qf-w-w ~ V-^c is obtained. 

A further piece of evidence for the alternating meson-glueball pattern is provided by looking 
at the excited pion spectrum, which we show in Fig. [3| The alternating pattern was unveiled 
by Glozman [27], suggesting that states degenerate with the pion might not be identified with 
glueballs. Remarkably, the states generating doublets with pion states are /o(600) o -7ro(140), 
/o(1370) o 7ro(1300), /o(1710) o 7ro(1800), /o(2100) o 7ro(2070), and /o(2330) o 7ro(2360), 
whereas the other scalar states /o(980), /o(1500), /o(2020) and /o(2200) are not degenerate with 
other mesons with light u and d quarks. Our analysis is reinforced by this observation. 

As a matter of fact, fitting the pion 7r(140) as the n = state of the Regge spectrum requires 
strong departure from a simple linear trajectory, = an + wlq. One may improve on this by 
using the holographic connection and a mixed boundary condition at z = determined by fixing 



the mass of the ground state mg using Eq. (17) together with Eq. (20) for the harmonic oscillator 
case, Eq. (19). This procedure ensures the orthogonality between all states and implements 
linearity for large n. This can be done for the ground states uiq = m7r, rra^, m/g. The fit to all 
states yields a = 1.37GeV^ and the mass spectra (in GeV) 

vr (Regge) (0.140, 1.260, 1.730, 2.092, 2.400, . . .) vr (PDG) (0.140, 1.300, 1.812, 2.070, 2.360) 
<T (Regge) (0.527, 1.297, 1.750, 2.106, 2.411, . . .) a (PDG) (0.600, 1.350, 1.724, 2.103, 2.321) 
/o(Regge) (0.977, 1.513, 1.906, 2.232, 2.517, . . .) /o(PDG) (0.980, 1.505, 1.992, 2.189) 

yielding l/-ya'i/'o(0)/V'o(0) = —3.1, —14.9, and 0.2, respectively. Note the large and small values 
for the a and /o cases, which suggests that these boundary conditions are very close to the 
Dirichlet and Neumann cases. Chiral symmetry breaking corresponds to the different tt and a 
values. 

7. Scalar dominance and heavy pions 

Hadronic matrix elements of the energy-momentum tensor, the so-called gravitational form 
factors (GFF) of the pion and nucleon, correspond to a dominance of scalar states in the large- A'^c 
picture, as {u{p) is a Dirac spinor) 

(vr(p')ieivr(p)) = (23) 
{Nip'mNip)) = u{P>ip)T.^-^^^^f^^ (24) 



where the sum rules grnnrfn = 1 28 and Mat = YlndnNNfn f29] hold. Unfortunately, the 
lattice QCD data for the pion |30| and nucleon (LHPC f311 and QCDSF |32| collaborations), 
picking the valence quark contribution, are too noisy as to pin down the coupling of the excited 
scalar-isoscalar states to the energy-momentum tensor. Nevertheless, useful information con- 
firming the (Regge) mass estimates for the c-meson can be extracted |16| through the use of 
the multiplicative QCD evolution of the GFF in the valence quark momentum fraction, {x)u+d, 
as seen in deep inelastic scattering or on the lattice at the scale (j, = 2 GeV. For the pion and 
nucleon GFFs we obtain the fits 

= 0.52(2), = 445(32) MeV, (25) 

= 0.447(14), m„ = 550;^80MeV. (26) 

Assuming a simple dependence of rrifj on m^r, 

m^(m^) = m^ + c(m^-m^^phys) , (27) 



8 



Ruiz-Arriola printed on October 14, 2011 



m^[MeV] 


[MeV] 


m„ [MeV] 


nia- [MeV] 


m„ [MeV] 




GFF 


GFF 


(gg)^-dynam. 


(gg) ^-quench. 


230 


580(190) 


620(100) 


— 


400(30) 


342 


630(190) 


660(90) 




720(20) 


478 


710(200) 


730(90) 




1000(20) 


318 


620(190) 


650(90) 


468(50) 




469 


700(190) 


720(90) 


936(13) 




526 


739(200) 


750(90) 


1066(13) 





Table 1. Scalar monopolc mass obtained from the nucleon gravitational form factors, extracted from the 
{qq) components obtained by LHPC 31 and QCDSF 32 and compared to the lattice calculation using 
the tetraquark {qq)'^ probing fields, both for the dynamical and quenched fermions 



33 



yields = 550l^^|jMeV and c 
on the choice of the lattice data 



31 



0.951^, 
or ^2^, 



or = 600l?goMeV and c = 0.8(2), depending 
respectively. Note that c is close to unity. Higher 
quark masses might possibly clarify whether or not the state evolves into a glueball or a meson. 
For a (qq)"' system one expects nio- — )• 2nmq + const at large current quark mass ruq. The data 
[32] are too noisy to see the difference, although for the largest pion masses we 
~ TTT-TT, as it simply corresponds to the gg-component of 0. We observe, however, 
~ 500MeV our results are not far away from the recent lattice calculation using 



from |31j or 
see that rria 
that for ttItt 
the tetraquark probes 



33 , which provide ma ~ 2m7r for the largest pion masses as they 



should (see Table [T]). From this viewpoint, and unless operator mixing is implemented, the 
nature of the state is predetermined by the probing operator. 



8. Dimension-2 condensates 



One of the problems of the large- A^c Regge models 35 and their holographic relatives 36 37 



is that they may contradict expectations from the OPE, as they involve dimension-2 condensates. 
For instance, the OPE for the Hqq^q'^) correlator in Eq.(12) gives corrections 0{q^), while the 
0{q'^) terms are missing [21] . This yields a one to one comparison: 



Co 



J n 



27r2 V a 



n 



2_2 



/3(a) 



a 



(28) 
(29) 

(30) 



Equation (28) requires infinitely many states, while Eq. (|29j) suggests a positive and non- 
vanishing gauge- invariant dimension-2 object, C2 = i J d'^xx'^{Q{x)@), which is generally non- 
local, as it should not appear in the OPE. Note that C2 > 0, hence is non-vanishing for a finite 
number of states. The infinite Regge spectrum of Eq. Q with Eq. (|28[ ) may be modeled with a 
constant /jq = = 0{Nc) whereas fn- = 0{y/Nc) goes as Eq. (14) and yields a convergent 
and hence positive contribution. Naively, we get C2 = 00. However, C2 may vanish, as required 
by standard OPE, when infinitely many states are considered after regularization. The use of 
the ("-function regularization [16^|34J gives 



^2 = 1™, E fnMsinf = E fn- + //o (1/2 



/a) 



(31) 



Then C2 = for ruf^ > ya/^ = 810(40) MeV, a reasonable value to 0{1/Nc) In any case, 
the important remark is that while at the OPE level C2 = vanishes for trivial reasons, at the 
Regge spectrum level some fine tuning must be at work. 
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